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ABSTRACT 

Context. The forthcoming new-generation radio telescope SKA (Square Kilometre Array) and its precursors will provide 
a rapidly growing number of polarized radio sources. 

Aims. Our analysis looks at what can be learned from these sources concerning the structure and evolution of magnetic 
fields of external galaxies. 

Methods. Recognition of magnetic structures is possible from Faraday rotation measures (RM) towards background 
sources behind galaxies or a continuous RM map obtained from the diffuse polarized emission from the galaxy itself. We 
constructed models for the ionized gas and magnetic field patterns of different azimuthal symmetries (axisymmetric, 
bisymmetric and quadrisymmetric spirals, and superpositions) plus a halo magnetic field. RM fluctuations with a 
Kolmogorov spectrum due to turbulent fields and/or fiuctuations in ionized gas density are superimposed. Assuming 
extrapolated number density counts of polarized sources, we generated a sample of RM values within the solid angle 
of the galaxy. Applying various templates, we derived the minimum number of background sources and the minimum 
quality of the observations. For a large number of sources, reconstruction of the field structure without precognition 
becomes possible. 

Results. Any large-scale regular component of the magnetic field can be clearly recognized from RM data with the help 
of the criterium. Under favorite conditions, a few dozen polarized sources are enough for a reliable result. A halo 
field with a vertical component does not affect the results of recognition. The required source number increases for 
small inclinations of the galaxy's disk and for larger RM turbulence. A flat number density distribution of the sources 
can be overcome by more sensitive observations. Application of the recognition method to the available RM data in 
the region around M 31 indicates that there are significant RM contributions intrinsic to the background sources or 
due to the foreground of the Milky Way. A reliable reconstruction of the field structure needs at least 20 RM values 
on a cut along the projected minor axis. 

Conclusions. Recognition or reconstruction of regular field structures from the RM data of polarized background sources 
is a powerful tool for future radio telescopes. Measuring RM at frequencies around 1 GHz with the SKA, simple field 
structures can be recognized in galaxies up to about 100 Mpc distance and will allow to test dynamo against primordial 
or other models of field origin. The low-frequency SKA array and low-frequency precursor telescopes like LOFAR may 
also have good RM sensitivity if background sources are still significantly polarized at low frequencies. 

Key words. Methods: statistical - techniques: polarimetric - galaxies: magnetic fields - galaxies: spiral - galaxies: 
individual: M 31 - radio continuum: galaxies 



1. Introduction 

Magnetic fields play an important role for the structure and 
evolution of galaxies. Radio synchrotron emission provides 
the best tools for measuring the strength and structure of 
galactic magnetic fields (Beck et al. 1996). However, the 
sensitivity of present-day radio observations only allows de- 
tailed studies for a couple of nearby galaxies (Beck 2005). 

Rotation measures ( RM) towards polarized background 
sources located behind spiral galaxies can trace regular 
magnetic fields in these galaxies out to large distances, even 
where the synchrotron emission of the galaxy itself is too 
weak to be detected (Han et al. 1998; Gaensler et al. 2005). 
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However, with the sensitivity of present-day radio tele- 
scopes, the number density of polarized background sources 
is only a few sources per solid angle of a square degree, so 
that only M 31 and the LMC, the two angularly largest 
galaxies in the sky, could be investigated so far (Han et al. 
1998; Gaensler et al. 2005). 

Future high-sensitivity radio facilities will be able to 
resolve the detailed magnetic field structure in galaxies 
by observing their polarized intensity and RM directly. 
Furthermore, they will be able to observe a huge number of 
faint radio sources, thus providing the high density back- 
ground of polarized point sources. This opens the possibility 
of mapping the Faraday rotation of polarized background 
sources towards nearby galaxies and study their magnetic 
fields. The detection of large-scale field patterns (modes) 
and their superpositions would strengthen dynamo theory 
for field amplification and ordering, while the lack of RM 
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patterns would indicate a primordial origin or fields struc- 
tured by gas flows (Beck 2006). Mapping the RM of the 
diffuse polarized emission of galaxies at high frequencies 
is restricted to the star-forming regions where cosmic-ray 
electrons are produced. At low frequencies the extent of 
synchrotron emission is larger due to the propagation of 
cosmic ray electrons, but Faraday depolarization increases 
much faster with wavelength, so that the diffuse polarized 
emission is weak. On the other hand, RM towards polarized 
background sources are generated by regular fields plus ion- 
ized gas, which extend to larger galactic radii, as indicated 
by the existing data sets (Han et al. 1998; Gaensler et al. 
2005), and are less affected by Faraday depolarization. 

A major step towards a better understanding of galactic 
magnetism will be achieved by the Square Kilometre Array 
(SKA, www.skatelescope.org) and its pathfinders, the Allen 
Telescope Array (ATA) in the US, the Low Frequency Array 
(LOFAR) in Europe, the McerKAT in South Africa, and 
the AustraUa SKA Pathfinder (ASKAP). The SKA will be 
a new-generation telescope with a square kilometre collect- 
ing area, a frequency range of 70 MHz to 25 GHz with con- 
tinuous frequency coverage, a bandwidth of at least 25%, 
a field of view of at least 1 square degree at 1.4 GHz, and 
angular resolution of better than 1 arcsecond at 1.4 GHz. 
The SKA is planned to consist of three separate arrays: a 
phased array for low frequencies (70-300 MHz), a phased 
array for medium frequencies (300 1000 MHz), and an ar- 
ray of single-dish antennas for high frequencies (1-25 GHz). 

"Cosmic magnetism" is one of six Key Science Projects 
for the SKA with the plan to measure a grid of more than 
10^ RM data towards polarized sources over the whole sky 
(Gaensler et al. 2004). This will allow measurement of the 
evolution of magnetic fields in galaxies from the most dis- 
tant to nearby galaxies and in the Milky Way. Large-scale 
field structures should be observable via RM mapping of 
the diffuse polarized emission in galaxies (Beck 2006). Even 
more promising are RM measurements towards polarized 
background sources, which are discussed in this paper. 

2. Models of Faraday rotation towards background 
sources in galaxies 

A model for RM towards polarized background sources 
needs the following ingredients: the distribution and sign 
of the regular and random field of the galaxy, the distribu- 
tion of electron density, the density and fluxes of polarized 
background sources, and their internal RM. The measured 
RM value towards each background source is the sum of all 
rotation measure contributions along the line of sight, in- 
cluding the Galactic foreground rotation and any intcirnal 
RM of the background source itself. The contribution of 
the Galactic foreground has to be subtracted from the ob- 
served RM values before any further analysis, but this con- 
tribution is presently known only for a coarse grid of back- 
ground sources on large scales of a few degrees (Simard- 
Normandin et al. 1981; Johnston-HoUitt et al. 2004). RM 
data of background sources near the plane of the Milky 
Way reveal fluctuations by several lOOradm"^ on scales 
of a few arcminutes (Brown et al. 2007) . The contribution 
of the Galactic RM can be accounted properly from the 
density of RM detected with the LOFAR or SKA in a 
slightly larger region surrounding each galaxy. The inter- 
nal RM of the background sources will average out if the 
number of background sources used for the analysis is large 



enough. Furthermore, a culling algorithm can be used to 
remove background sources with high internal RM values. 
Hennessy et al. (1989) and Johnston-Hollitt et al. (2004), 

for example, removed sources with measured RM exceeding 
n-times the standard deviation from the mean (or median) 
RM. 

2.1. Models for the regular magnetic field in the disk 

The mean-field a-Q dynamo model is based on differential 
rotation and the a-effect (Beck et al. 1996). Although the 
physics of dynamo action still faces theoretical problems 
(e.g. Brandenburg & Subramanian 2005), the dynamo is the 
only known mechanism able to generate large-scale coherent 
(regular) magnetic flelds of spiral shape. These coherent 
fields can be represented as a superposition (spectrum) of 
modes with different azimuthal and vertical symmetries. In 
a smooth, axisymmetric gas disk the strongest mode is the 
one with the azimuthal mode number m = ( axisymmetric 
spiral field), followed by the weaker m = 1 {bisymmetric 
spiral field), etc. (Elstner et al. 1992). These modes cause 
typical variations in Faraday rotation along the azimuthal 
direction in the galaxy disk (Krause 1990). In flat, uniform 
disks the axisymmetric mode with even vertical symmetry 
(SO mode) is excited most easily (Baryshnikova et al. 1987), 
while the odd symmetry (AO mode) dominates in spherical 
objects. The timescale for building up a coherent field from 
a turbulent one is ~ 10^ yr (Beck et al. 1994). 

Most galaxies reveal spiral patterns in their polarization 
vectors, even flocculent or irregular galaxies. The field ob- 
served in polarization can be anisotropic or regular. RM is a 
signature of regular fields, while anisotropic fields generate 
polarized emission, but no Faraday rotation. Large-scale 
RM patterns observed in several galaxies (Krause 1990; 
Beck 2005) show that at least some fraction of the mag- 
netic field in galaxies is regular. The classical case is the 
strongly dominating axisymmetric field in the Andromeda 
galaxy M 31 (Berkhuijsen et al. 2003; Fletcher et al. 2004). 
A few more cases of dominating axisymmetric fields are 
known (e.g. the LMC, Gaensler et al. 2005), while dom- 
inating bisymmetric fields are rare (Krause et al. 1989). 
The two magnetic arms in NGC 6946 (Beck & Hoernes 
1996), with the field directed towards the galaxy's center 
in both, arc a signature of superposed m = and m = 2 
modes. However, for many of the nearby galaxies for which 
multi-frequency observations are available, angular resolu- 
tions and/or signal-to-noise ratios are still too low to re- 
veal dominating magnetic modes or their superpositions. 
The regular field of the Milky Way reveals several reversals 
in the plane, but the global structure is still imclear (Han 
et al. 1997; Brown et al. 2007, Sun et al. in press). The 
regular field near the Sun is symmetric with respect to the 
Galactic plane. 

For this paper, we restrict our analysis to the three 
lowest azimuthal modes of the toroidal regular magnetic 
field: axisymmetric spiral (ASS, m = 0), bisymmetric spiral 
(BSS, m = 1), quadrisymmetric spiral (QSS, m = 2), and 
the superpositions ASS+BSS and ASS+QSS. All modes are 
assumed to be symmetric with respect to the disk plane (S- 
modes) . The maximum regular field strength is assumed to 
be 5 fxG for all modes, consistent with typical values from 
observations (Beck 2005). 

The radial extent of regular magnetic field is not known 
yet. In NGC 6946, its scalelength is at least 16 kpc (Beck 
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2007). RM measurements of polarized sources behind M 31 
indicate that the regular field strength out to at least 15 kpc 
radius is similar to that in the inner disk (Han et al. 1998). 

2.2. Models for the magnetic field in the halo 

Galactic dynamos operating in thin disks also generate 
poloidal fields that extend far into the halo. They arc about 
one order of magnitude weaker than the toroidal fields in 
the disk (Ruzmaikin et al. 1988) and hence are irrelevant 
for our models. Polarization observations of nearby galax- 
ies seen edge-on generally show a magnetic field parallel to 
the disk near the disk plane, but recent high-sensitivity ob- 
servations of several edge-on galaxies like M104, NGC 891, 
NGC 253, and NGC 5775 show vertical field components 
that increase with increasing height z above and below the 
galactic plane and also with increasing radius, so-called X- 
shaped magnetic fields (Heesen et al. 2005, 2007; Soida 
2005; Krause et al. 2006; Krause 2007). Such a vertical field 
can be due to a galactic wind that may transport the disk 
field into the halo. Dynamo models including a wind out- 
flow show field structures quite similar to the observed ones 
(Brandenburg et al. 1993). 

The detailed analysis of the highly inclined galaxy 
NGC 253 allowed a separation of the observed field into an 
ASS disk field and a vertical field (Heesen et al. 2007). The 
tilt angle between the vertical field and the disk is about 
45°. This value corresponds approximately to the tilt angles 
at large z observed in other edge-on galaxies. 

Assuming energy-density equipartition between mag- 
netic fields and cosmic rays, the scale height of the total 
magnetic field is about 4 times larger than the scale height 
of synchrotron emission of typically about 1.8 kpc (Krause 
2004). As the degree of linear polarization increases with 
height above the disk midplane, the scale height of the reg- 
ular field may be even larger. 

2.3. Thermal electron density models 

In this paper we model the ionized gas of the disk by a 
Gaussian distribution in the radial and vertical (z) direc- 
tions. A scale height of Ikpc is adopted from the models 
of the free electron distribution in the thick disk of the 
Milky Way fitted to pulsar dispersion measures (Gomez 
et al. 2001; Cordes & Lazio 2002). The electron density 
near the galaxy center is assumed to be 0.03 cm~^, the 
same as in the Milky Way models. 

The radial scalelength of the density Ue of free electrons 
in the Milky Way of «15 kpc (Gomez et al. 2001) is much 
larger than the scalelength of radio thermal emission in 
nearby galaxies, e.g. about 4 kpc in NGC 6946 (Walsh et al. 
2002), 3 kpc in M 33 (Tabatabaei et al. 2007), and 5 kpc 
for the Ha disk of the edge-on galaxy NGC 253 (Heesen, 
priv. comm.). As the thermal emission scales with n^, the 
scalelength of Ue is twice larger, but still smaller than what 
is quoted for the Milky Way. As a compromise, we use a 
Gaussian scalelength of 10 kpc in our model. 

The distribution of ionized gas follows that of the opti- 
cal spiral arms. The arm-interarm-contrast of the Ha emis- 
sion is between 2 and 10 in typical spiral galaxies like 
NGC 6946 (see the smoothed Ha map shown in Frick 
et al. 2001). On the other hand, regular magnetic fields 
observed in polarization are strongest in interarm regions 
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Fig. 1. The two coordinate systems used in this paper: I, b 

in the sky plane and s along the line of sight, and x, y, z in 
the model galaxy plane and perpendicular to it. The dashed 
lines indicate the thickness of the disk. 



(Beck & Hoernes 1996; Frick et al. 2001) with a similar arm- 
interarm-contrast. As Faraday rotation in NGC 6946 is also 
highest in the interarm regions (Beck 2007), the azimuthal 
variation in Faraday rotation, determined by the product of 
regular field strength and electron density, is mostly deter- 
mined by the azimuthal variation of the regular field. We 
conclude that the azimuthal variation in electron density 
can be neglected in our models. 

In summary, we adopt a Gaussian dependence of rig on 
galactic radius r (in the galactic plane) and on height z 
above the disk midplane: 

ne(r,2) = rioexpj- (^^^ jexpj-^^^ |, (1) 

where no=0.03 cm~^ is the electron density near the galaxy 
center, ro=10 kpc is the radial scalelength of the Ug distri- 
bution in the galaxy model, and h=l kpc is its vertical 
scalelength. 

2.4. Faraday rotation models 

Two Cartesian coordinate systems are used: the coordinate 
system l,b,s in the sky plane with the same orientation as 
the galaxy (/, b are the distances counted from the center of 
the observed galaxy along the major and minor axes of the 
projected disk, and s corresponds to the direction of the 
line of sight) and the system x,y, z attached to the galaxy 
(x, y are in the galactic midplane and z is perpendicular to 
it). The coordinates x and I of both systems coincide and 
i is the angle between coordinates z and s (or y and b, see 
Fig. 1), called the inclination angle of the galactic plane 
(where i = 0° means facc-on). 

The galactic magnetic field is modeled as a superposi- 
tion of a regular part with a simple azimuthal symmetry 
Bm and a random part Bturbi which describes the contri- 
bution of large-scale galactic turbulence: 

B{r) = B„{r) + Bturbir), (2) 
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Fig. 2. Modeled strength of the regular magnetic field multiplied by the thermal electron density in the galactic midplane 
(i.e. in the plane x,y) for the field modes m = (ASS), m = 1 (BSS), and to = 2 (QSS) for B = 1. 



where r is the radius vector in galactic coordinates, which 
will be written in Cartesian coordinates (x,y,z) or cylindri- 
cal ones (r = \/ x"^ +2/^: 4' = arctany/x, z). 

First, the regular part of magnetic field is supposed to be 
purely horizontal {Bz = 0), parameterized by the pitch an- 
gle p and the intensity Bm of the corresponding azimuthal 
mode TO (we only consider to = 0, 1, 2, using the notation 
ASS, BSS, and QSS, see Fig. 2). The BSS mode has two 
reversals along azimuthal angle, the QSS mode has four re- 
versals. All modes are symmetric with respect to the disk 
plane. The strength of the magnetic field in cylindrical co- 
ordinates is defined as 

Bm{r.,(j},p) = Scos ^TO (y^^ - + </'o^^ tanh (3) 

where B is the field amplitude (strength) and the az- 
imuthal phase of the mode. The tanh term is introduced 
to suppress the field near the center of the galaxy where 
the spiral field is strongly twisted. Equation (3) does not 
include any decrease in the field strength with radius r and 
height z, because the decay in regular field strength is much 
slower than that of electron density (see Eq.(l)). 

The RM is proportional to the product of the density 
of thermal electrons rie and the component of the regular 
magnetic field parallel to the line of sight. Note that 
we cannot separate the contributions of B and rig using the 
RM data alone and that we can reconstruct only the prod- 
uct ( B\\ Tie). As we neglect any azimuthal (e.g. spiral) struc- 
ture in TLf,, the small-scale structure in RM is determined 
by B, whereas the general radial and vertical decrease in 
RM is determined by ng. The distribution of this field in 
the galactic midplane for to = 0, 1, 2 (ASS, BSS, and QSS) 
is shown in Fig. 2. 

Second, we consider a model in which a vertical field B^ 
exists and increases with height z, according to the mag- 
netic field structure observed in the halos of nearby edge-on 
galaxies (see Sect. 2.2). We suppose that the modulus of the 
magnetic field follows again Eq. (3), but the field lines are 
tilted with respect to the plane. The tilt angle x is defined 
as 

X = Xo tanh tanh — , (4) 

where xo = t^/^ is the limit of the tilt angle in our model 
achieved at r w 7'o/3, z k, 2h. The orientations of magnetic 
field lines for the ASS model are shown in Fig. 3. 




Fig. 3. Orientation of the halo magnetic field. The box 

sketches the region of ionized gas r < rQ, —h < z < h. 



The magnetic field component i^y parallel to the line 
of sight includes the turbulent (random) magnetic field as 
well as the regular one, so that 

/oo 
{Brn{r,p) + Bturb{r))n,{r)ds. (5) 
-oo 

The turbulent part of magnetic field Bturb should describe 
the three-dimensional random field with given spectral 
properties in the whole range of scales. To avoid full 3-D 
simulations of the random vector field and its contribution 
to RM, we prefer to model this contribution of the irregu- 
lar part of magnetic field directly in the RM maps. Thus 
we write 

/oo 
B,^{r,p)n,{r)ds+miturb{lM, (6) 
-oo 

where we add the random part of the galactic field projected 
to the sky plane, i.e. integrated along the line of sight. 

WAturb{lib) is produced in a way that allows for the 
spectral and spatial distribution of galactic turbulent fields. 
Vogt & EnBlin (2005) show that by observing a turbulent 
magnetic field with a three-dimensional Fourier power-law 
spectrum like |B(fc)p ~ fc" one obtains an RM map, the 
two-dimensional spectrum, from which another power law 
follows I R"M(fc_L)p ~ fc_L^ with the slope (3 ^ a. Here the 
hat is used for the notation of the Fourier transform of 
any function /(fe) = j f{x)e^''^dx, k is the 3-D wave vec- 
tor and k± is the 2-D wave vector, defined in the plane of 
the RM map. It is supposed that the magnetic spectrum 
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Fig. 4. Modeled RM maps (in radm~^) for an inclination angle of i = 10° generated for a pure BSS field pattern, for 
B = 5 iiG and no=0-03 cm~'^ (left), a random turbulent field KMturb for AKMturb = 30 radm~^ (middle), and their 
superposition (right). 



is Kolmogorov- like {EB{k) ~ k'^\B{k)\^ ~ k'^/^, hence 
a = —11/3) at large wave numbers (on scales smaller than 
the thickness of the galactic disk) and is growing at small 
wave numbers as Esik) ~ k^ (then a = 0). Thus we gen- 
erate a random map in the Fourier space with the required 
spectrum of the form: 

Here $o defines the amplitude of the turbulent part of the 
generated map, and fcg is the wave number corresponding to 
the turbulence macroscale (fco = 5 is adopted here to reach 
a realistic spectrum with a maximum on scales correspond- 
ing to the typical size of supernova remnants). Hereafter we 
return to physical space and obtain the map RMt,irb(^&) 
as the product of the random map <f>(Z, h) and the galactic 
"electron thickness" : 

KMturb{l.h)^'^{lMDM{l,b). (8) 

Namely, DM {I, b) is the dispersion measure for the given 
line of sight: 

/oo 
ne{r)ds. (9) 

The amplitude of the signal in this map is normalized by the 
parameter $0j so that the r.m.s. ARM(„r6 of the random 
field KyVturb adopts some given value. We will consider 
three cases: weak turbulence (ARMt„rb = 10 radm~^), 
moderate turbulence (A RMt„rb = 30 radm~^), and strong 
turbulence [l^'RMturb = 50 radm~^). Note that, due to the 
normalization, A RMt„r& does not depend on the amplitude 
of DM and hence does not increase with increasing inclina- 
tion, although such a behavior may occur in real galaxies 
(see Sect. 6). 

Figure 4 shows an example of the RM map calculated 
for a galactic field observed almost face-on (inclination an- 
gle i = 10°), simulated for a pure BSS field pattern (left), 
the map of turbulent field RMt,jr6, and the resulting RM 
map for the same BSS model superimposed by a moderate 
level of turbulence with a "realistic" spectrum (right). 



3. Number density of polarized background sources 

Future high-sensitivity radio facilities (LOFAR, SKA) will 
be able to observe faint radio sources, thus providing a high- 
density background of polarized sources. This opens a new 
possibility of mapping the Faraday rotation of polarized 
background sources towards nearby galaxies and studying 
their magnetic fields. In this section we estimate the den- 
sity of the polarized sources towards a typical spiral galaxy 
placed at different distances from the observer and study 
the errors associated with rotation measures and their de- 
pendence on the telescope specifications. 

The differential number counts of polarized sources has 
been recently studied using the deep polarimetric observa- 
tions at 1.4 GHz of the European Large Area ISO Survey 
North 1 region (ELAIS Nl) as a part of the DRAG Planck 
Deep Fields project (Taylor et al. 2007). They show that 
the observed distribution of number counts in differential 
representation (i.e. numbers per flux density interval) is al- 
most flat down to polarized flux densities of « 0.5 mJy. 
They also show that the extrapolated differential number 
counts of polarized sources decreases down to the 0.01 mJy 
level. The corresponding cumulative source number counts 
is shown in Fig. 5, normalized to the solid angle of 1 deg^, 
which roughly corresponds to the area covered by the inner 
part of the nearest galaxy M 31 (Han et al. 1998). 

Gbservations with the SKA (assuming 10^ m^ collect- 
ing area, 50 K receiver system temperature and ^vjv — 
AA/A = 0.25 relative bandwidth) will achieve an r.m.s. 
noise of 0.1 /xJy/beam area and 0.01 ^Jy/beam area at 
1.4 GHz (A=21 cm) within one hour and 100 hour integra- 
tion time, respectively (Garilli & Rawlings 2004). 

We extrapolate the source counts at A=21 cm from Pg = 
0.5 mJy to the limiting flux density Pmin of 0.05 /xJy, as- 
suming a power-law relation between the cumulative num- 
ber counts per deg^ (A^n) and polarized flux density (P): 

Nn{>P) = f^^[^ \ (10) 

where = 16 deg~^ is the cumulative source number at 
P = Pg = 0.5 mJy. The exponent 7 is expected in the 
range between 0.7 and 1.1 ("pessimistic" and "optimistic" 
number counts). The number count per solid angle of 2ro x 
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Fig. 5. The dotted line sliows the cumulative count N of 
polarized sources brighter than the polarized flux density P 
at A=21 cm, taken from Taylor ct al. (2007), extrapolated 
to 0.01 mJy polarized flux density and normalized to source 
numbers per 1 deg^ solid angle on sky. The possible range 
of extrapolated number counts per galaxy solid angles of 
1 deg^ are shown as thick solid lines, 0.013 dcg^ as dashed 
lines, and 0.00013 deg^ as dot-dashed lines. These solid an- 
gles correspond to distances of a galaxy (with ro = 10 kpc) 
at 1.15 Mpc, lOMpc, and 100 Mpc. The two different lines 
at each distance refer to slopes 7 of the extrapolated source 
numbers of 0.7 (lower line) and 1.1 (upper line). The two 
thin vertical lines represent the 5ap detection limit for 
sources observed at A=21 cm with the SKA within 100 
hours (left) and 1 hour (right) observation time. 



2ro around a typical face-on spiral galaxy with a radius of 
ro at a distance D is 



N = Non„ ( ^ 



(11) 



where fig ~ 



{2ro/Df{180/TTf is the solid angle of the 
galaxy, given in square degrees. 

In Fig. 5 we show the predicted number counts per 
typical spiral galaxy (rg = 10 kpc, 20 x 20kpc^ area) ob- 
served at a distance of 1.15 Mpc (1 deg^ solid angle), 10 Mpc 
(0.013dcg2) and 100 Mpc (0.00013deg2). With the "quick- 
look" sensitivity of 0.01 mJy, the SKA will already detect 
within a few minutes about 1000 polarized sources towards 
the nearest spiral galaxy M 31. With the best achievable 
SKA sensitivity of P « 0.05 /iJy, about 5x10'' polarized 
sources behind M 31 and hundreds to tens of polarized 
background sources will be detected behind a typical galaxy 
at a distance between 10 and 100 Mpc. Only for the pes- 
simistic case (7 = 0.7), the source munber becomes very 
small for galaxies near 100 Mpc distance. 

The instrumental Icr error in rotation measure A RM of 
each simulated radio source of polarized flux density (P) in 
the cumulative number counts is 



ARM 



2A2^' 



A the center wavelength, and AA/A the relative bandwidth 
of the observations. If the instrumental r.m.s. noise in polar- 
ized flux density is almost constant across the whole band, 
(T;^ is determined by the errors of polarization angles in the 
frequency channels at the smallest and largest wavelength, 
which for simplicity are assumed to be the same. Then, 



= V2^, 
2P' 



(13) 



where dp is the instrumental r.m.s. noise in polarized flux 
density per frequency channel and P is the polarized flux 
density of the source. Hence, 



ARM = 



where 



P' 



A = 



2V2A2 AA- 



(14) 



(15) 



The "sensitivity parameter" A only depends on instrumen- 
tal specifications. It becomes smaller for lower noise ap 

(larger integration time of the telescope), larger relative 
bandwidth, and longer observation wavelength. For the ex- 
pected performance of the SKA at A=21 cm (see above), 
A « 3 /iJy rad m^^ /^/t where t is the observation time 
in hours. For a more general discussion of RM errors, see 
Appendix A of Brentjcns & de Bruyn (2005). 

The maximum error of RM is reached at the limiting 
flux density Pmin of the survey. For P^ 



50-p one derives 



A RM^Qa,(P^j^) 



14A2 AA ' 



(16) 



which is independent of noise ap and of observation time. 
For the SKA ARM^ax is « 6 rad m"^. For comparison, 
the VLA (using two bands centered at A=22 cm and 18 cm, 
see Sect. 4.5) gives AKMmax ~ 10 rad m"^. 

Typical maximum values of RM arc RM « 40 rad m"^ 
for a galaxy inclined by i = 10° (almost face-on) and 
RM w 200 rad m"^ for i = 90° inchnation (edge-on). 
With the SKA centered at A=21 cm, the maximum error 
A KMmax is significantly smaller than the typical maximum 
values of RM for edge-on and face-on nearby galaxies. 

From Eqs. (14) and (10-11) one may recover 



ARM = gArV7^ 



where 

Q = A{Nong) 



-1/7 p-i 
-'0 • 



(17) 



(18) 



The "observation parameter" Q only depends on obser- 
vational characteristics. It defines ARM for the bright- 
est source statistically expected within the solid angle of a 
given galaxy by a given telescope (i.e. for N=l in Eq. (17)). 
It varies with noise Up and with distance D as OpD"^!^ . 
The value of Q lies in the range between « 0.001 rad m"^ 
and « 10 rad mT^ for Op between 0.01 /iJy and 0.03 mJy 
and the range of distances from 1 Mpc to 100 Mpc (Fig. 6). 
Prom Eqs. (11) and (18), the parameter Q directly relates 
to the maximum number of sources N^ax (above the 5crp 



(12) detection limit Pmin) per sohd angle of the galaxy: 



where as is the la error (in rad) of the total difference in N^„^(P^ 
polarization angles </> within the observed wavelength range. 



= ( 10x/2A2 



AA 



Q- 



(19) 
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Fig. 6. Relation between instrumental noise and distance 
to the galaxy (see Eq. (18)) for different values of the "ob- 
servation parameter" Q — 0.001,0.01,0.1,1,10 rad m^^ 
(from left to right) drawn for pessimistic and optimistic 
slopes, 7 — 0.7 (dashed Hne) and 7 = 1.1 (solid Unc) for 
A — 21cm, 6X/X — 0.25. The dotted line corresponds to 
the observed slope 7^1 (see Fig. 5). The numbers of 
sources per solid angle of the galaxy corresponding to the 
above Q values are N^ax — 462, 92, 18, 4, 1 for 7 = 0.7 and 
Nmax = 15411, 1224, 97, 8, 1 for 7 = 1.1 (Eq. (19)). 



4. Recognition of regular galactic field structures 
from a limited sample of background sources 

In this section we describe and use the minimization 
method to test the magnetic field models in the disk and 
halo of a galaxy. Application of this method to recognize 
the magnetic field structure in M31 is also presented. 

4.1. Method 

Suppose that N polarized background sources are observed 
with an "observation parameter" Q towards a galaxy with 
radius rg at distance D. Given the Faraday rotation values 
of the model (RM„iod) caused by the regular field (RM^eg) 
and the turbulent field in the galaxy, we want to recognize 
the structure of the regular magnetic field, such as ASS, 
BSS, QSS, or their combinations, and determine the best- 
fit parameters of the model (amplitude, spiral pitch angle, 
and azimuthal phase) and their standard deviations. 

For a fixed inclination angle of the galaxy, we simulate a 
"template" of Faraday rotation values KMreg,n{^i, ■■■,£,d) 
of N points (with coordinates (/„, bn)) for the regular field 
only where are the free parameters of the magnetic field 
model: the amplitude and spiral pitch angle of the mag- 
netic field (number of free parameters c? = 2 for ASS), the 
azimuthal phase (d = 3 for BSS and QSS), and the rel- 
ative amplitude q for the ASS field in combination with 
the BSS or QSS field where the BSS or QSS field has the 
amplitude of 1 — g (cZ = 4 for the combination). For the 
same set of points we simulate the "observed" KMrnod. n 
(Eq. (5)) which takes both the regular and turbulent mag- 
netic fields into account. By varying the free parameters we 
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Fig. 7. Factor of effective points vs inclination angle i due 
to geometry of the disk of thermal electrons 



minimize the normalized value between the "observed" 
and "template" RM maps. 



1 



N 

E 



(RM, 



RM, 



(ARM 



2 

noise,n 



ARM 



.(20) 



and estimate the best-fit parameters for each model. 
(A KMnoise and A HMturb are the instrumental noise 
(Eq. (17)) and the noise of the random field, respectively.) 
The coefficients Wn define the weight of each point Wn — 
DM{ln, bn)/DM{0, 0), hence the contribution of each point 
to the minimization procedure. The weight of the point Wn 
is defined by the dispersion measure DM(/„,6„) (Eq. (9)) 
along the line of the sight. Introducing this weight is needed 
to balance the contribution of strong RM sources with 
small error, which are located at the galactic periphery, 
where the low thermal electron density does not allow to ob- 
tain reliable information about the galactic magnetic field. 
The total value of Wn gives the weighted (effective) num- 
ber of points Nm = Wn- This number is different from 
the real number of points N taken per 2ro x 2ro area. It 
also depends on the inclination angle i (see Fig. 7). This 
dependence is similar to cosi, but it is not equal to 1 at 
i = 0° because of the Gaussian distribution of the ther- 
mal electrons and not equal to at i = 90° due to the 
finite thickness of the disk. The inclination angle 45° yields 
iVu, = O.AAN. This means that in order to obtain at least 
Nyj = 4 we need a minimum of A'^ = 10. 

This method is suited for simple field structures and a 
relatively small number of RM points allowing the limited 
number of parameters to be fitted. In the case of compli- 
cated field structures and/or a large number of RM points, 
a full reconstruction should be attempted (Sect. 5). 

4.2. Fitting the simulated data 

Throughout this paper we use typical values of = lOkpc, 
h ^ Ikpc, p = 20°, ric = 0.03 cm^^ in Eq. (1) and 
B = 5 fiG in Eq. (3). We start with recognition of an ASS 
field as the simplest model. For a fixed i = 45°, we sim- 
ulate the RM values, noise distribution ARM, and tur- 
bulent components of N uniformly distributed background 
sources. We determine the best-fit parameters p' and B' 
by applying the minimization technique (Eq. (20)). 
Repeating this procedure for several thousand realizations 
we calculate the mean pitch angle {p'), amplitude (i?'), and 
their standard deviations {6p' and 6B'). The dependence 
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Fig. 9. Required asymptotic values N* and Q* to recognize 
an ASS-type field. The points correspond to 5p' = 3° for 
three different levels of turbulence (from bottom to top in 
each triplet: fO radm~^, 30 radm"^ and 50 radm^-^) and 
for different inclination angles (15°, 45°, 75°) for fixed h = 
1 kpc, p = 20°. Filled symbols correspond to the optimistic 
case (7 — 1.1), open symbols to the pessimistic one (7 — 
0.7). 



Fig. 8. ASS field model: mean values (figures in left col- 
umn) and standard deviations (figures in right column) of 
pitch angles p' and amplitudes B' (given in fiG) for fixed 
values oi h = 1 kpc, p = 20°, B = 5 fiG, i = 45° and assum- 
ing the pessimistic case of number counts (7 = 0.7). N is 
the real number of points observed within the solid angle of 
the galaxy and Q is the observation parameter (Eq. (18)). 



of the statistical properties of the fitted parameters on the 
number N of sources per solid angle of the galaxy and on 
Q are shown in Fig. 8. 

The simulated value of p' is in good agreement with the 
value p = 20° used for the model calculation over a wide 
range of N and Q. However, the standard deviation of 6p' is 
quite small only for low values of Q and for large numbers 
N. The amplitude B' is well defined, too. The behaviour of 
SB' is the same as for 5p' . The standard deviation diagrams 
show that there is a certain value of Q below which better 
observations do not improve the result of fitting. The same 
is valid beyond a certain value of N. For example, for the 
level of Sp' — 3° no more than 36 sources observed with 
Q = 0.2 are needed, or no better values than Q ^ 0.05 
when observing 21 sources. Hence, the fitting problem can 
be characterized by two asymptotic values N* and Q* , pro- 
vided that the accuracy of the fitting is fixed. (For the given 
example one obtains N* = 21 and Q* = 0.2.) These two 
parameters define the required quality of observations. 

In Fig. 9 we show the dependence of asymptotic param- 
eters N* and Q* on the level of turbulence and pitch angle 
of the observed galaxy for optimistic and pessimistic values 
of 7. The required level of accuracy is fixed at Sp' = 3°. 

We conclude that, in the case of the pessimistic slope, 
one needs observations with smaller Q, hence smaller noise 
(Tp, which requires a longer integration time. An increase in 
the number of sources (for a fixed noise value Cp) does not 
help. This can be explained by the fact that for low values 
of 7 additional sources are fainter, have higher errors, and 
do not improve the result. The optimistic evaluation of 7 
essentially increases the values of Q*, but slightly affects 



the values of N* . As expected, the result of fitting strongly 
depends on the level of turbulence of the galactic field: in 
the case of intense turbulence one needs more points and 
smaller Q (longer observations or a more sensitive instru- 
ment). The dependence on turbulence becomes dramatic 
for weakly inclined (face-on) galaxies - the reliable fitting 
requires a huge number of sources under any value of 7. 

The same simulations are performed for the BSS and 
QSS models. An additional free parameter (the phase 0) 
is accounted for, but, in general, the results are similar to 
those obtained for the ASS model (see Figs. 10 and 11). 
Comparing Figs. 9 and 11 (or Figs. 8 and 10), we con- 
clude that the bisymmetric field can be easier recognized 
(it requires fewer points and higher Q for the same accu- 
racy of pitch angle definition) and weakly depends on the 
turbulent component of magnetic field. The BSS field has 
better chances of being recognized in slightly inclined galax- 
ies. However, with values of , the BSS model has a higher 
probability of being false than the ASS model (see Fig. 12). 

The dependence of the fitting accuracy on the incli- 
nation angle is shown in Fig. 13. We learn that face-on 
galaxies always pose more problems for recognition of the 
field structure. The accuracy of the ASS or BSS recogni- 
tion increases monotonically with inclination angle due to 
stronger line-of-sight components of the regular field. For 
inclination angles > 70° , the number of effective points de- 
creases. Furthermore, RM becomes sensitive to fluctuations 
of the regular field, e.g. due to spiral arms (see Sect. 6), so 
that the reconstruction method (Sect. 5) is recommended 
for strongly inclined galaxies. 

To illustrate the potential of identification of different 
models of the regular magnetic field, using a set of tem- 
plates, we present in Table 1 the results of recognition un- 
der a moderate level of galactic turbulence {AKMturb — 
30 rad m~^) for a sample of 35 sources. The models of mag- 
netic field structures with fixed parameters of p = 20°, 
00 = 0°, and B = 5 fiG are used. The assumed level of 
galactic turbulence is AKMturb = 30 rad m"^, the incli- 
nation of the galaxy's disk is 45°, the slope of the source 
counts 7 = 0.7 ("pessimistic case"), and the "observation 
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Fig. 10. BSS field model: Mean values (figures in left col- 
umn) and standard deviations (figures in right column) of 
fitted pitch angles p', amplitudes i?', and azimuthal phases 
(fo for fixed /i = 1 kpc, p = 20°, (/iq = 0°, and i = 45°, and 
assuming the pessimistic case of number counts (7 = 0.7). 

is the real number of points observed within the solid 
angle of the galaxy and Q the observation parameter. 



parameter" is Q = 0.04. The statistics for each case has 
been done over 220 random distributions of 35 points within 
the solid angle of the galaxy and 30 different turbulent 
fields for each distribution (in total 220x30=6600 realiza- 
tions). We adopt a significance level of 98% for the rejection 
of a model, corresponding to the probability function 
Cl{v/2,x^v/2) < 0.02 where = 7V„ - d is the weighted 
number of degrees of freedom (Eq. (20)) and N^^ ~ 27 for 
i — 45° (Fig. 7). A model is rejected at the 98% confidence 
hmit if > 1.7 (see Table 1). 

In all cases the model is clearly identified as the tem- 
plate with the lowest value of i almost reaching the limit 
of 1.0 for a perfect fit. The pitch angle is reproduced almost 
precisely, with an uncertainty that decreases from ±2° for 
the ASS field to ±1° for the QSS field. The accuracy of the 
field amplitude is about 10% for all field modes. 
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Fig. 11. Required asymptotical values N* and Q* to rec- 
ognize a BSS-type field. The points correspond to Sp' — 3° 
for three different levels of turbulence (from bottom to top 
in each triplet: 10 radm^^, 30 radm^^ and 50 radm~^) 
and different inclination angles (15°, 45°, 75°) for fixed 
h = Ikpc, p = 20°. Filled symbols correspond to opti- 
mistic case (7 = 1.1), open symbols to the pessimistic one 
(7 = 0.7). 
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Fig. 12. Probabilities calculated from foi' ASS (left) and 
BSS (right). 



4.3. Influence of vertical magnetic field 

We also tested models including a vertical magnetic field 
in the halo, as described in Sect. 2.4 and shown in Fig. 3. 
We find that such a vertical magnetic field in the halo only 
slightly changes the RM map even at small inclination an- 
gles. 

The infiuence of the vertical component of magnetic 
field on the recognition procedure is studied by the fol- 
lowing test. We calculate the RM map for a model with a 
vertical magnetic field and try to recognize it by a template 
having only horizontal magnetic fields. As shown in Fig. 14, 
the fitted pitch angles and amplitudes of the magnetic field 
are in the same range as without a vertical magnetic field 
for galaxies with inclination angles i > 15°, but large devi- 
ations occur for smaller inclinations. The reason becomes 
clear from Fig. 3: For small inclinations the line-of-sight 
components of the vertical field are larger than those of the 
field in the plane and hence distort the results, while at 
larger inclinations the vertical field becomes unimportant. 
We conclude that the regular magnetic field structure of 
spiral galaxies (at least for inclinations of i > 15°) can be 
recognized well by a magnetic field template without a ver- 



10 



Stepanov et al.: Reconstruction of magnetic field structures 



Table 1. The "x^ Championship" 



template 
model J, 


ASS 


BSS 


QSS 


P' 


B' 




p' 




B' 


x' 


p> 




B' 


x' 


ASS 


20±2 


5.0±0.3 


1.1±0.4 


7±18 


2±96 


4.6±1.8 


21±6 


23±48 


-6±41 


5.2±0.8 


13±3 


BSS 


-22±47 


0.6±0.7 


13±3 


21±2 


14±32 


5.1±0.4 


l.litO.3 


8±14 


10±52 


4.0±1.9 


11±3 


QSS 


-25±39 


0.7±0.6 


12±3 


8±23 


-8±95 


3.5±1.3 


9±3 


20±1 


-3±16 


4.9±0.5 


1.3±0.3 


ASS+BSS 


19±6 


2.5±0.4 


3.8±1 


19±16 


13±107 


3.2±1.1 


6±2 


24±47 


2±37 


2.9±0.7 


6±2 


ASS+QSS 


18±6 


2.6±0.4 


3.4±1 


4±14 


16±95 


3.0±1.5 


8±3 


10±35 


10±33 


3.3±1.2 


6±2 



template — > 
model J, 


ASS+BSS 


ASS+QSS 


P' 


g' 




B' 


x' 


P' 


g' 




B' 


x' 


ASS 


20±2 


0.9+0.1 


-3+98 


5.6+0.3 


1.2+0.4 


20+3 


0.9+0.1 


7+52 


5.7+0.4 


1.2+0.3 


BSS 


21±2 


0.0+0.1 


9+26 


5.0+0.4 


1.2+0.4 


8+14 


0.1+0.1 


2+50 


5+2 


10+3 


QSS 


2±19 


0.1+0.1 


-21+95 


3.7+1.2 


10+4 


20+1 


0.0+0.1 


0+6 


5.1+0.3 


1.2+0.3 


ASS+BSS 


22+3 


0.5+0.1 


18+46 


5.1+0.4 


1.2+0.4 


18+10 


0.6+0.1 


6+50 


4.0+0.8 


3.7+0.9 


ASS+QSS 


15+10 


0.7+0.1 


-3+96 


4.1+0.8 


3.3+1.3 


21+2 


0.5+0.1 


6+27 


4.9+0.5 


1.2+0.3 



Parameters of fitted pitch angles p' and field amplitudes B' , azimuthal phases 0oi s^nd values of the best-fit template (first row). 
The models of field structures are given in the first column. For the superposition templates (second part of the Table), q' is the 



relative amplitudes of the ASS mode, with q 
rejected are marked in bold. 



: 0.5 assumed in the model. The models with x < 1-7 which cannot be statistically 
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Fig. 13. Accuracy (standard deviation) of the fitted pitch 
angle (upper panel) and the magnetic field amplitude (lower 
panel) for the ASS model (open circles, black lines) and 
the BSS model (triangles, red fines) vs inclination angle i. 
Weak turbulence - dashed lines, strong turbulence - solid 
lines. Fixed parameters of the model are the pitch angle 
p = 20° and the field amplitude B ^ 5 /iG. Each template 
is based on 35 RM points. 
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Fig. 14. Fitted pitch angle (upper panel) and the magnetic 
field amplitude (lower panel) for the ASS model (open cir- 
cles, black lines) and the BSS model (triangles, red lines), 
both including a vertical field, recognized with templates 
without a vertical field, plotted vs inclination angle i. Weak 
turbulence - dashed lines, strong turbulence - solid lines. 
Fixed parameters of the model are the pitch angle p = 20° 
and the field amplitude B = 5 fiG. Each template is based 
on 35 RM points. 



tical field, even if the galaxy has such a vertical magnetic 
field. 

4.4. Effect of position uncertainty 

Our method assumes that the position of the galaxy cen- 
ter known is precisely. However, this may not be the case 



for distant galaxies for which only low-resolution optical 
images from surveys are available. A 2" uncertainty corre- 
sponds to 0.5 kpc at 50 Mpc and 1 kpc at 100 Mpc distance. 

To test the effect of position uncertainties, we shifted 
the templates by 0.5 kpc and 1 kpc with respect to the 
models of type ASS, BSS, and QSS. The best fits have the 
same parameters of the field structures as in Table 1, but 
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Fig. 15. RM source sample of M31 (from Han et al. 
(1998)). Numbers give RM (in radm^^), after subtrac- 
tion of the mean RM from the Galactic foreground. Isolines 
show the weight function used for ^ . Coordinates are given 
in kpc. 



the values for the ASS field increase to 1.3 and 1.8 for 
0.5 kpc and 1 kpc shifts, respectively, to 1.5 and 2.4 for the 
BSS field, and to 2.2 and 4.3 for the QSS field. This tells us 
that the disturbing effect of a position uncertainty increases 
with mode number. A typical 2" uncertainty is uncritical 
for ASS-type fields even at 100 Mpc distance, while QSS- 
type fields can no longer be recognized at distances greater 
than about 50 Mpc. 

4.5. Application of tlie recognition method to M 31 

Han et al. (1998) measured 22 RM values from polarized 
background sources within a solid angle of about 1 deg^ in 
the inner part of M 31 (Fig. 15). Their VLA observations 
at A=22 cm and at A=18 cm (AA/A « 0.015) reached a 
sensitivity limit of 0.3 mJy/beam area. The source number 
agrees well with the source counts from Fig. 5. Han et al. 
(1998) claimed that their RM values are consistent with an 
ASS field in M 31, but did not include a statistical test. 

Observations of the diffuse polarized emission and its 
RM from M 31 has been described by an ASS field with 
8° — 19° pitch angle, 4 amplitude (assuming ric — 
0.03 cm~'^ and 1 kpc pathlength), with a weaker QSS field 
superimposed in the inner part of the galaxy (Fletcher et al. 
2004) , while earlier results indicated a superimposed weaker 
BSS field (Sofue & Beck 1987). 

Firstly, we estimate how representative the sample of 
Han et al. (1998) is. We use the same models with fixed 
parameters and assumptions as described in Sect. 4.2 and 
repeat the same test as described above (see Table 1) for 
different model fields, but using for the recognition only the 
coordinates of the 22 sources (shown in Fig. 15) and the 
observed instrumental errors A IUs/L„oise for each source, 
but we used the RM values obtained from the field models. 
The results in Table 2 show that the sample of 22 sources 
is large enough to clearly recognize the field structure with 
simple symmetry (if it exists). Models with > 2.7 are 
rejected at the 98% significance level, which corresponds to 
a weighted number of degrees of freedom of = N^, — d'^b, 
where ~ 8 (for the inclination of M 31 of z = 75°) and 
the number of free parameters is d = 2, 3, or 4. 



N 


temp. 


P' 


q' 


<l>'o 


B' 






ASS 


-8±31 


- 


- 


0.4+0.2 


5.3+1.2 




BSS 


12±16 


- 


46+121 


2.2+2.6 


5.0+1.1 


22 


QSS 


3±45 


- 


-10+53 


1.7+2.1 


5.0+1.2 




A+B 


6±16 


0.1+0.1 


-51+115 


2.7+1.9 


5.7+1.4 




A+Q 


44+34 


0.3+0.1 


-39+38 


2.6+1.7 


5.1+1.2 




ASS 


-16+22 






0.4+0.2 


4.6+1.2 






1 -t-91 




oomyz 


o. / mo.D 


o.omi.o 


20 


QSS 


9+36 




-12+52 


2.4+2.7 


5.1+1.2 




A+B 


-17+17 


0.4+0.2 


34+89 


2.+1.5 


6.3+1.7 




A+Q 


8+32 


0.2+0.2 


-28+50 


2.8+2.1 


6.4+1.6 




ASS 


-3+17 






0.5+0.2 


2.8+0.9 




BSS 


8+15 




17+120 


2.8+2.4 


3.3+1.2 


20 


QSS 


28+39 




-12+42 


2.6+3.1 


3.3+1.1 




A+B 


3+19 


0.3+0.2 


9+114 


3.1+2.3 


3.8+1.4 




A+Q 


19+34 


0.3+0.1 


-21+49 


3.6+2.1 


3.7+1.2 



Parameters of fitted pitch angles p', field amplitudes B' , az- 
imuthal phases t/iQ, amplitude ratios q (in case of a superposi- 
tion), and of the best fit, using 22 and 20 RM values observed 
in the M 31 regime by Han et al. (1998). For 22 and 20 values 
(8 and 6.5 effective points), models with > 2.7 and x^ > 3.1, 
respectively, are rejected at the 98% significance level. 



Secondly, we use the observed data to apply our recog- 
nition model. The analysis was performed in a first run for 
all 22 points. In a second run we exclude the two sources 
with the highest RM deviations from the template model 
field (those with +57 rad m~^ and -17 rad m~^), follow- 
ing the culling algorithm as described in Sect. 2 in order 
to reduce the influence of strong internal RM of the back- 
ground sources. As shown in Table 3 none of the two runs 
give a clear result. Using all 22 sources, all templates give a 
similar value of x^. Taking only 20 sources (middle part of 
Table 3), a small preference for the ASS model is indicated. 
However, all models are rejected at the 98% significance 
level (x^ > 2.7). As the simple field structures could be 
clearly identified in our test using the same set of sources 
(Table 2), the results of Table 3 seem to indicate that the 
structure of the regular magnetic field of M 31 is more com- 
plicated than suggested by our templates. However, the 
analysis of the diffuse emission by (Fletcher et al. 2004) 
is statistically reliable and does not require higher modes 
with significant amplitudes. We conclude that the measured 
RM values towards the background sources contain a signif- 
icant, intrinsic contribution or a fluctuating contribution of 
the Milky Way foreground, which can be subtracted only if 
a larger number of background sources becomes available. 

Thirdly, to account for the RM intrinsic to the sources 
or from the Galactic foreground, we increase all RM errors 
by an arbitrary factor of three. The x^ values of the best 
fit decrease and the ASS model becomes slightly preferred 
statistically (see bottom part of Table 3) , in agreement with 
the results from diffuse RM. 



5. Magnetic field reconstruction from RM grids of 
inclined galaxies 

In the previous section the possibility of recognizing the 
large-scale structure of the regular magnetic field using a 
small number of sources has been discussed. In this one 
we switch to the case of a large number of sources, which 
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Table 2. The "x^ Championship" for M31. 



template 
model I 


ASS 


BSS 


QSS 


p> 


B' 


x' 


p> 




B' 




p> 




B' 




ASS 


20±1 


5.0±0.2 


1.3±0.6 


14±1 


128±4 


8.8±0.5 


84±7 


62±6 


-5±3 


4.2±0.2 


71±10 


BSS 


0±58 


O.OitO.l 


76±8 


20±1 


0±3 


5.0±0.2 


1.4±0.6 


15±2 


11±35 


9.1±0.9 


34±5 


QSS 


-2±31 


0.4±0.2 


14±3 


19±11 


-122±77 


3.7±3.5 


9±2 


20±1 


0±3 


5.1±0.5 


1.5±0.5 


ASS+BSS 


30±5 


1.5±0.2 


13±2 


14±1 


165±4 


4.8±0.5 


20±3 


58±8 


-10±31 


1.9±0.4 


13±2 


ASS+QSS 


20±3 


2.6±0.2 


4.8±1.4 


14±1 


123±17 


4.7±0.4 


25±4 


58±19 


-6±10 


2.5±1.4 


24±4 



template — > 
model J, 


ASS+BSS 


ASS+QSS 


p' 


q' 




B' 


x' 


p' 


<?' 




B' 


X 


ASS 


20+1 


0.9+0.1 


15+112 


5.6+0.3 


1.7+0.8 


20+1 


0.9+0.1 


5+58 


5.6+0.3 


1.8+0.7 


BSS 


20+1 


0.0+0.1 


0+5 


5.0+0.3 


1.9+0.9 


15+2 


0+0.1 


7+39 


9+1 


43+5 


QSS 


14+13 


0.0+0.1 


-121+82 


5.2+3.9 


11.4+2.5 


20+1 


0.0+0.1 


0+1 


5.1+0.6 


1.5+0.5 


ASS+BSS 


20+2 


0.5+0.1 


5+27 


5.0+0.3 


1.5+0.7 


25+5 


0.4+0.1 


0+12 


4.4+1.5 


8+2 


ASS+QSS 


20+4 


0.7+0.1 


83+125 


3.7+1.2 


4.2+1.7 


20+1 


0.5+0.1 


2+13 


5.1+0.6 


1.5+0.5 



Parameters of fitted pitch angles p' and field amplitudes B' , azimuthal phases 0q, and of the best fit for different model fields 
(first column), recognized by a template (first row) based on the existing sample of 22 sources in the M 31 area. Models with 
> 2.7 (marked in bold) are rejected at the 98% significance level. 



allows us to reconstruct the structure of the galactic field, 
without precognition of a given set of simple models. The 
method discussed below is especially powerful for resolving 
the detailed structure of the magnetic field of a strongly 
inclined galaxy, provided an RM grid with strong enough 
resolution is available. 

Observations reveal a grid of RM(/,6) ± ARM(;,6) 
where the standard deviation ARM(^,6) characterizes the 
observational noise. Note again that here, in contrast to the 
recognition method discussed in the previous section, we do 
not separate RM(Z, b) into its regular and turbulent parts, 
but instead merely reconstruct the existing field structure, 
trying to suppress the infiuence of random signals (instru- 
mental or generated in the galaxy). The Faraday rotation 
measure is proportional to the product (iJyne). Ignoring 
the problem of the and separation, we consider the 
problem of reconstructing the magnetic field ( rie) from 

the observed RM grid. We suppose that ( B|| rie) can be 
written in the form: 



Bn n,. 



f{x,y) ri{z) 



This allows us to describe the decrease in |RM| in the z- 
direction by one term ri{z) and to state the problem as the 
reconstruction of the 2-D scalar function f{x,y) (see Fig. 1 
for the coordinate system). 

The aim is to reconstruct the field f{x,y) from the RM 
grid. The rotation measure grid, even in the case of noiseless 
measurements, is a convolution of the required field f{x,y) 
with some function g{y), defined by the vertical structure 
of the galactic field. For the case = 0, the function g(y) 
is defined only by the function r](z). Hence, 



RM(/, b) = 

= C / f f{x, y) ri{z) 6 (z~ ^ ] dydz = 



— oo ./ — oo 
oo 



smi 



= C" / fix,y)g{b/smi~y)dy, 
where 



(22) 



9{y) 




Fig. 16. The function f{x, y) (left) and the RM grid (right) 
for a BSS galaxy model (pitch angle p — 20°) inclined at 
i = 70° . 



where i is the inclination angle and C a constant. 

The standard solution for a problem like Eq. (22) is 
based on the Fourier deconvolution. Actually, 



(21) f{x,k)^ 



RM(Z,fc) 

m 



(23) 



77 (ycosz) dz, 



where k is the wave number. The hat denotes the Fourier 
transform of the corresponding function. The Fourier de- 
convolution gives the exact solution of Eq. (22) , but requires 
some regularization procedure (Heinz & Martin 1996) to 
avoid the amplification of high spatial frequency noise, es- 
pecially if the spectrum of the kernel function g decreases 
rapidly with k. The often-suggested Gaussian form of the 
galactic profile rj{z) leads to a kernel that is "bad" in this 
sense. We do not discuss the regularization problem here 
in detail. The simplest algorithms use some kind of high- 
frequency cut-off in the Fourier space, while more sophis- 
ticated ones introduce subtle filtration, based for example 
on wavelet algorithms (Johnstone et al. 2004). 

The main problems of using Eq. (23) arise from the 
limited resolution of the RM grid (small number of points) , 
noise, and unknown kernel function g (namely, one does not 
know the vertical distribution function rj{z)). 

As an example let us consider a strongly inclined galaxy 
with a pure BSS-type magnetic field. The required function 
/(x, y) for this case and the simulated RM grid for the in- 
clination angle i = 70° are shown in Fig. 16. The problem 
can be illustrated by reconstructing the one-dimensional 
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Fig. 17. The profiles f{y) along or parallel to the minor 
axis of the projected galaxy disk in the sky plane and the 
corresponding profiles of RM(6) for a BSS-type galaxy, in- 
clined at i = 70°: a; = (thick curve), a; = 3 kpc (dashed), 
and X = 7.5 kpc (thin). 



profiles f{y) for given cross-sections x = const. Figure 17 
shows the profiles of f{y) for three different values of x and 
the corresponding three profiles of RM(6) at the inclina- 
tion i = 70°. To use the deconvolution (23) one needs the 
function r]{z). Even if a Gaussian shape is adopted, one free 
parameter (the half-thickness of the ionized gas disk h) re- 
mains. We denote by h' the supposed value of h used for 
the deconvolution. 

Figure 18 shows the profiles f{y) reconstructed from 
RM(&) by deconvolution with different assumed values for 
the half thickness h' . The simulation is done for a set 
of points corresponding to 20 sources across the galaxy 
(which seems to be quite realistic in the frame of the SKA 
project). This figure illustrates that the result of recon- 
struction strongly depends on the choice of h' . An under- 
estimated value of the thickness of the ionized disk leads 
to insufficient reinforcement of small-scale structures - the 
deconvolution merely stretches the RM projection. In con- 
trast, an overestimation of h' implies an overrated smooth- 
ing and results in a disproportional amplification of high 
spatial frequencies (small scales). According to the last 
panel of Fig. 18, the oscillations of f{y) essentially exceed 
the real solid angle of the galaxy. The correct value h' = h 
allows restoration of two positive maxima in the central 
part of the minor axis profile a; = 0, which were completely 
hidden in the RM projection. This example tells us that a 
slight underestimate of h' is more acceptable than an over- 
estimate. Figure 18 also shows that the central profile at 
a: = is the most informative one for the reconstruction, 
while the outer profiles are more similar. 
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Fig. 18. The profiles f{y) reconstructed using three differ- 
ent values of the supposed galactic thickness: h' — 0.75 kpc 
(top), h' ^ 1 kpc (middle), h' = 1.25 kpc (bottom). The 
correct value is ft. = 1 kpc. The profiles are given for a; = 
(thick curve), a; = 3 kpc (dashed), and x = 7.5 kpc (thin). 



6. Discussion 

Up to now the analysis of the magnetic field structure in ex- 
ternal galaxies from RM observations is still in a rudimen- 
tary state. There are only two nearby galaxies behind which 
a limited sample of polarized background sources were ob- 
served (M 31 and LMC). The impressive projects of large 
radio telescope arrays (LOFAR and SKA) promote the in- 
vestigations of future possibilities. We considered two kinds 
of analyzes that may become possible: the detailed recon- 
struction of the structure of the regular magnetic field in 
a number of nearby galaxies from a dense grid of RM val- 
ues from background sources in the galaxy field, and the 
recognition of the main symmetry of the large-scale regular 
magnetic field from the RM of a number of background 
sources. We studied this problem based on current esti- 
mates of the increasing number density of polarized sources 
with improving telescope sensitivity. We took an optimistic 
and a pessimistic extrapolation of the existing cumulative 
number counts into account. 
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In the first part of the paper, we studied the possibility 
of recognizing the large-scale regular component of the mag- 
netic field with a simple symmetry using a straightforward 
testing by templates, which correspond to axisymmctric, 
bisymmctric, and quadrisymmctric modes, or their super- 
positions, and applying the critcrium for evaluating the 
reliabihty of the different modes. We showed that, if a sym- 
metric part indeed exists in the analyzed field (even mixed 
with a relatively strong turbulent field), dozens of sources 
provide a good chance of a reliable recognition. These tem- 
plates can be even successfully applied if the disk field is ac- 
companied by an X-shaped vertical field as found in nearby 
galaxies observed almost edge-on. The result strongly de- 
pends on the quality of the observations, the flux density 
of the sources, and their number density distribution. 

Our results show that regular magnetic fields of sim- 
ple azimuthal symmetry (modes), like axisymmetric (ASS), 
bisymmctric (BSS), and quadrisymmctric spiral fields 
(QSS), and superpositions of these, can be recognized with 
RM values from a number of background sources depend- 
ing on the inclination of the galaxy's disk, the level of RM 
turbulence in the galaxy, and the slope of the number den- 
sity distribution of the polarized sources. The fitted pa- 
rameters are the spiral pitch angle, the amplitude of the 
regular field strength, and the phase (for non-axisymmetric 
fields). The typical accuracies achieved for the fitted pa- 
rameters are ±(1° — 4°) for the pitch angle and ±10% for 
the field amplitude. Under the favorite conditions, 10 — 30 
sources are already reliable enough to recognize the field 
structure. However, the unknown internal RM of the back- 
ground sources themselves may require a larger number of 
measured sources in order to statistically average out the 
contribution of internal RM, or to be able to apply a culling 
algorithm, as mentioned in Sect. 2. 

The regular field of a higher order of symmetry can be 
easier recognized; i.e. it requires fewer points and shorter 
observation time for the same accuracy of the fitted pitch 
angle definition and is less affected by the turbulent compo- 
nent of RM. The BSS and QSS fields seem to have better 
chances of being recognized in slightly inclined galaxies. 
However, the smaller required N* and the larger number of 
parameters for the higher mode decreases the reliability of 
the recognition. 

The dependence on turbulence becomes dramatic for 
weakly inclined (almost face-on) galaxies a reliable fitting 
requires a huge number of sources. This problem may occur 
especially for small, slowly rotating galaxies with a high 
rate of star formation, where the large-scale regular field is 
expected to be weaker than the turbulent field (Beck et al. 
1996). Interactions or ram pressure may also enhance the 
turbulent field and mask any large-scale pattern in the RM 
distribution. 

In our models we did not consider deviations of the reg- 
ular field structure from simple azimuthal modes, e.g. due 
to spiral arms or additional field reversals. These would 
lead to further RM fluctuations that may become large for 
strongly inclined galaxies and hence decrease the accuracy 
of the flts. Strong RM fluctuations have been observed near 
the plane of the Milky Way (Brown et al. 2007). 

For our test case M 31, the best case of a dominat- 
ing field mode (ASS) so far, the available RM data from 
polarized background sources confirm the results from the 
analysis of RM data from the diffuse emission of M 31 
(Sect. 4.5) if the RM errors are increased to account for 



the RM contributions intrinsic to the sources or from the 
foreground in the Milky Way. A larger number of sources 
is needed for proper subtraction of these effects. However, 
it is possible that higher field modes exist and will prevent 
a clear result even in case of a greater source number. 

The detection of the ASS mode in M 31 from RM of 
background sources also indicates that this mode is sym- 
metric with respect to the plane (mode SO) because RM 
in a strongly inclined galaxy hosting an AO mode should 
reverse above and below the plane. The RM amplitude of 
polarized sources behind a SO-type field should be twice 
larger than that from the RM map of the diffuse polarized 
emission from the galaxy itself, which has to be confirmed 
with more sensitive data. 

For many galaxies no strongly dominating field mode 
can be expected. A high value, which does not decrease 
with increasing source number, would indicate that the field 
structure is more complicated than a superposition of a few 
simple modes or that the turbulent field is stronger than the 
regular field. In this case, our reconstruction method should 
be applied, because it does not need a "precognition" tem- 
plate but needs a higher density of RM sources and hence 
deeper observations. 

Observations provide the 2-D projection of the product 
(B|| rie). Interpretation of this projection can by impeded 
by a vertical magnetic field , but the contribution of this 
field is not crucial in the case of symmetrically tilted fields. 
If the inclination angle increases, the RM grid is entirely 
dominated by the horizontal field, but another problem 
arises: the smoothing of detailed structures (of the size of a 
spiral arm) due to their superposition in a highly inclined 
projection. In other words, the function f(x,y) describing 
the field structure in the galaxy's plane is convolved by a 
smoothing function, defined by the vertical profile of the 
galaxy. Provided a universal shape of the profile exists, one 
can apply the deconvolution technique to restore the struc- 
ture of the magnetic field in the galaxy's plane. In principle, 
this method works for any field structure. 

The reconstruction method is superior for strongly in- 
clined galaxies (about 70° and more) and is succcssfiil if 
a large sample of sources is available. At least about 20 
sources are required for one cut along the minor galac- 
tic axis (or ^ 1200 sources within the solid angle covered 
by the galaxy) to reconstruct the spiral arms in a spi- 
ral galaxy. The number of points increases if the scale of 
reconstructed details decreases, following the well-known 
Kotclnikov-Nyquist relation. 

Note that the reconstruction procedure can be applied 
to the RM obtained from the diffusive polarized emission of 
the galaxy itself, but this would require the deconvolution 
of a more complicated integral equation than Eq. (22). On 
the other hand, the recognition method also works for a 
continuous RM map obtained from the diffuse polarized 
emission. 

Finally, we wish to point out that little is known about 
the statistical properties of the RM contributions intrin- 
sic to the background sources and of the contribution of 
small-scale RM inhomogcneities in the foreground of the 
Milky Way. Large, nearby galaxies will suffer most from 
fluctuations in the foreground. Another uncertainty are de- 
polarization effects within the background sources, or in 
the foreground of the galaxy or the Milky Way, if the an- 
gular extent of a source is greater than the angular tur- 
bulence scale in the foreground medium. Depolarization of 
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background sources in the LMC were detected by Gaensler 
et al. (2005). The latter depolarization effect increases with 
distance because larger turbulence scales can contribute. 
Both depolarization effects, internal and external, depend 
on angular resolution and on observation frequency. No sta- 
tistical data are available yet. Further investigations in this 
direction may modify the number of sources required for a 
reliable recognition. 

7. Application to LOFAR and SKA 

The reconstruction of magnetic field structures of strongly 
inclined spiral galaxies is possible for a sample of > 1200 
RM soiirces. This would require a sensitivity of the SKA 
at 1.4 GHz of » 0.5 — 5 /xJy (or integration time less than 
one hour) for galaxies at distances of about one Mpc. The 
field structures of galaxies at about 10 Mpc distance can be 
reconstructed with tens to a hundred hours of integration 
time. 

The results for the recognition method presented in this 
paper are very promising for future observations with the 
SKA of background sources behind galaxies that are too 
faint to be observed directly by their extended polarized 
emission. Simple field structures can be recognized with the 
SKA at 1.4 GHz (21 cm) with an "observation parameter" 
Q K 1. This needs a o-p « 0.2 ^Jy for galaxies at 30 Mpc 
distance and (jp w 0.015 /iJy for galaxies at 100 Mpc dis- 
tance. These sensitivities can be achieved within 15 min and 
100 h observation time, respectively, assuming that polar- 
ization calibration of the SKA is possible down to such low 
flux-density levels. 

Let us estimate the number of spiral galaxies for 
which one can reconstruct or recognize the magnetic fields. 
Marinoni et al. (1999) estimate the mean number of galax- 
ies (~ 0.032 Mpc~^; the multi-attractor model was used 
to predict the distances to the galaxies) for a magnitude- 
limited all-sky optical sample of nearby ^ 5300 galaxies 
with recession velocities of cz < 5500 km s^^. The mean 
density of 2240 Sbc-Sd spirals is about 0.013 Mpc"^, which 
counts for about 60 spirals within a distance of 10 Mpc and 
« 60000 spirals within 100 Mpc. 

The recognition method raises requirements for the dy- 
namic range of the telescope for polarization measurements 
(polarization purity), which can be directly taken from 
Fig. 5. The detection of about 100 polarized sources in the 
solid angle of a galaxy at about 1 Mpc distance needs a dy- 
namic range of about 20 dB, while at 10 Mpc distance 30 dB 
may already be needed for a pessimistic value of 7, and even 
40 dB at 100 Mpc distance. Hence, reliable estimates of 7 
are essential to compute the required polarization purity. 

Observing at longer distances D requires increasing the 
observation time T according to T oc D^/t to obtain the 
same number N of sources per solid angle of the galaxy 
and the same observation parameter Q. Better knowledge 
of the slope 7 of cumulative number counts in the flux range 
accessible to the SKA is crucial for observations of distant 
galaxies. 

The existence of simple field modes in a majority of 
galaxies would give strong evidence of dynamo action. 
Failure to detect simple field structures would suggest that 
large-scale dynamos are unimportant in galaxies or that 
other processes like shearing or compressing gas flows de- 
form the field lines (Beck 2006). Complicated field struc- 
tures require application of the reconstruction method. 



which will become possible for a large number of relatively 
nearby galaxies with the SKA, hereby offering the chance 
to get good statistics about the magnetic field modes and 
their superpositions that may exist in nearby galaxies. 

Field recognition may turn out to be also useful at 
low frequencies. Equations (16) and (18) illustrate the ad- 
vantage of radio telescopes operating at long wavelengths. 
For example, LOFAR will provide a maximum bandwidth 
of 32 MHz in its lowband (30 80 MHz) and highband 
(110-240 MHz). With the large number of spectral chan- 
nels available for LOFAR, bandwidth depolarization will 
be negligible. The channel width determines the maxi- 
mum detectable RM. More severe is the maximum scale 
in Faraday depth beyond which an extended source with 
internal Faraday rotation (i.e. emission and rotation occur 
in the same volume) becomes depolarized (Brentjens & de 
Bruyn 2005); it varies with A^^„, is below 1 rad m~^ in the 
LOFAR bands and will limit low-frequency RM measure- 
ments of extended sources. 

LOFAR is planned to reach an r.m.s. noise (within 
one hour of observation time and with 4 MHz bandwidth) 
of 1 2 m.Jy/beam area in the lowband and better than 
0.1 mjy/beam area in the highband. This gives a sensi- 
tivity parameter for RM measurements of (Eq. (15)) of 
A PS 30 — 50 /xJy rad m~^ in both bands, compared to 
« 3 AtJy rad m'^ with the SKA at 1 GHz. LOFAR will al- 
low RM to be measured from sources at the detection limit 
with errors of A RM^aa; ~ 0.005 rad in the lowband 
and ARMtoox ~ 0.1 rad m~^ in the highband. However, 
the RM errors will be much greater under realistic condi- 
tions due to the limitations of correcting for ionospheric 
Faraday rotation. 

The low-frequency array of the SKA should reach 
an r.m.s. noise of 1.4 /xJy/beam area within one hour 
at 200 MHz (CariUi & Rawhngs 2004), hence A « 
1 /xJy rad m^^, much better than with the SKA at higher 
frequencies and better than with LOFAR, while the max- 
imum RM error of AKMmax » 0.1 rad m~^ is similar to 
that using the LOFAR highband. 

However, the number density of polarized sources at 
low frequencies may be much lower than at high frequencies 
that would increase the observation parameter Q (Eq. (18)) 
and hence reduce the number of RM points per solid an- 
gle of the galaxy. First, the overall degree of polarization is 
most probably lower due to Faraday depolarization effects. 
Second, the typical degree of polarization may decrease for 
more distant sources. Star formation was stronger in young 
spiral galaxies, and lobes of radio galaxies were smaller. 
Third, internal and external depolarization may signifi- 
cantly reduce the number of polarized sources (Sect. 6). 
For planning RM surveys at frequencies of 300 MHz and 
lower, much more needs to be known about the statistics 
of polarized sources at low frequencies. 

Another problem at low frequencies is thermal absorp- 
tion that may reduce the radio fluxes at frequencies below 
about 100 MHz for strongly inclined galaxies with high den- 
sities of the ionized gas. 

8. Conclusions 

Recognition of simple structures of regular magnetic flelds 
in galaxies, characterized by spiral pitch angle, amplitude 
of regular held strength, and azimuthal phase, can be reli- 
ably performed from a limited sample of RM measurements 
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towards polarized background sources. Applying templates 
to the field models, we could show that single and mixed 
modes can be clearly recognized based on the criterium. 
The quality of RM observations can be characterized by the 
number of detected polarized sources within the solid an- 
gle of the galaxy, varying with distance D and instrumental 
noise in polarization ap a.s N oa a~'^D~^, and the "obser- 
vation parameter" Q, which varies as Q oc UpD^/^ and also 
depends on wavelength and bandwidth of the observations. 
Field recognition is successful beyond a minimum source 
number N* and below a maximum observation parameter 
Q*. 

The slope 7 of the cumulative number density counts 
of polarized sources with flux densities P {N oc P~^) 
is actually known beyond P = 0.5 mJy from available 
source counts at 1.4 GHz. Extrapolation to smaller P yields 
0.7 < 7 < 1.1. 7 still needs to be determined, also at 
lower and higher frequencies. For the optimistic case of 
7 = 1.1, significant inclination of the galaxy's disk (i > 25°) 
and weak RM turbulence (due to turbulent magnetic fields 
and/or fluctuations of the ionized gas density), a few dozen 
polarized sources are already sufficient to reliably recognize 
a simple field structure of type ASS, BSS, or QSS. Higher 
modes (BSS and QSS) arc easier to recognize; i.e., they 
need k^ss soiirccs or yield smaller errors of the pitch angle 
for the same source number compared to the basic ASS 
mode. The requirement for a telescope's dynamic range in 
polarization (polarization purity) may exceed 30 dB in case 
of a fiat slope 7 of the cumulative number density counts 
of polarized sources. 

Vertical fields do not affect the results of the reconstruc- 
tion for galaxies with inclinations i > 15°. The method 
should not be applied for galaxies at smaller inclinations 
where RM from the disk field is small and a vertical field 
may dominate RM. 

Strong turbulence or small inclination need more po- 
larized sources for a statistically reliable recognition. In 
the case of small inclination angles, RM turbulence can 
get comparable to the RM amplitude of the regular field 
and the required number of polarized sources may become 
unreasonably large. As the required number of polarized 
sources per solid angle also increases for large inclination 
angles due to geometry, the recognition method works best 
for intermediate inclination angles. 

Uncertainties in the position of the galaxy's center of a 
few arcseconds arc even acceptable for galaxies at 100 Mpc 
distance if the field pattern is type ASS. For BSS-type fields 
accuracies of better than 2" are needed beyond 50 Mpc 
distance, and better than 1" beyond 50 Mpc distance for 
QSS-type fields. 

The application of the recognition method to the avail- 
able RM values within the solid angle of the galaxy M 31 is 
consistent with the analysis of the diffuse polarized emission 
if significant RM contributions intrinsic to the background 
sources or fiuctuating RM in the Milky Way foreground 
are assumed. The regular field of M 31 is type SO, i.e. ax- 
isymmetric spiral in the plane and symmetric with respect 
to the plane. 

Reconstruction of the field structure without precogni- 
tion needs a large number of sources. A reliable reconstruc- 
tion of the field structure is possible for at least 20 RM 
values on a cut along the projected minor axis. 



Future RM observations with the SKA at about 1 GHz 
will extend the possible targets for field recognition to 
about 100 Mpc distance, allow detailed statistical studies 
of the frequency of field structures, dependencies on galaxy 
properties, and test dynamo against primordial models or 
other models of field generation. About 60 spiral galaxies 
can be observed within a distance of 10 Mpc and w 60000 
spirals within 100 Mpc. 

Radio telescopes operating at low frequencies (LOFAR, 
ASKAP, and the low-frequency SKA array) may also be 
useful instruments for field recognition or reconstruction 
with the help of RM, if background sources are still signif- 
icantly polarized at low frequencies. 
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